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Abstract 

We study the spectral problems for the spatially periodic flows of inviscid incompressible 
fluid. The basic flows under consideration are the shear flows whose profiles oscillate on high 
frequencies. For such flows, we present asymptotic expansions of the unstable eigenvalues in the 
case when the limit spectral problem has multiple eigenvalues. 

Introduction 

The Euler equations of inviscid incompressible fluid admit simple steady solutions known as the 
shear flows. Such solutions (in the Cartesian coordinates) are written in the form 

V = (f/(i/),0,0) P = const 

where V denotes the velocity field, P denote the pressure, and U is referred to as the flow profile. 
Generally, function U can be chosen arbitrarily but we restrict our considerations within the case 
of smooth 2-77— periodic profiles. For such flows, we study 2D spatially periodic perturbations using 
the linear approximation. Let G denote the stream function for the velocity perturbation; hence 
the instant velocity of the perturbed flow is 

V + {Gy{x,y,t),-G^{x,y,t),0). 

Then the linear equation for perturbations takes the form 

f) fin 

AGt + U{y)^AG - U" — = 0, A = 9^ + (1) 

Since Eq. ([1]) is invariant with respect to the the translations both in t and in x, it is natural to 
seek for 

G{x,y,t) = e'^*e^°^$(y). 
This substitution leads us to the periodic spectral problem for the Rayleigh equation 

{a + iaU)(^-^-a'^^^-iaU"^ = 0; ^{y) = ^{y + 2tt). (2) 
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If there exists an eigenvalue a with Reo" > then there exists spatially periodic perturbation that 
grows exponentially in t and the basic shear flow is treated as unstable one. Following to [T] and 
[2], we look for such instability in the case of rapidly oscillating basic flow. 

Let us consider profile U{y) = W{s), where W - smooth 27r-periodic in s, and s = my with 
integer m (m >> 1). Let us seek for the solutions of ([2]) in the form 

m = e'^'^fis). (3) 

where n is integer and / is unknown 27r-periodic function. Substituting ([3]) into ([2]) yields equation 
for /: 

i/((T,e)/ = 0, (4) 

where operator H is defined by 

Hf = {a + iaW)f" - iaW" f+ 
+2ine{(T + iaW)f' - e^ic? + r?){a + iaW)f, 

accent means differentiation by variable s. We construct the asymptotic for the unstable eigenvalues 
a in the case of e = ^ — > 0. Our approach follows to that of ^ and In more details, let us 
write the problem Q in the form 

H^j + eFi/ + E^H^f = 0, (5) 

Hof = 0/" - 0'7; Hs) = CT + iaWis); (6) 

Hif = 2in(t>f'; ; (7) 

H2f = -{a^ + n^)(Pf;. (8) 

Assuming that a — > do when e — > we get 

^00/ = {^0 + iaW{s))f" - iaW'f = 0. (9) 

Eq. ([9j) has periodic solution 

f{s) = C{ao + iaW{s)) 

for every do G C. To make this solution unique (up to the constant factor C) one have to require 

7(ao) =< (do + iaWis))'^ >^ 0, (10) 

where <> denotes the averaging over the period. The analysis of discovers the unstable eigen- 
values under assumption (jlOp . (In other words, the eigenvalue of the limit problem is required to 
be simple.) In particular, this condition is satisfied for W{s) = sins provided RecJo ^ 0. Then do 
have to be selected using the solvability condition of the equation of the next approximation to the 
eigenf unction. For a generic profile, however, 7 has zeroes outside the real axis. For such ctq Eq. ([9]) 
has two independent solutions, i.e. ctq represents a multiple eigenvalue. In this paper, we develop 
the asymptotic in the case of multiple ctq. In particular, we prove that each zero of 7 generically 
gives rise to the unique branch of simple eigenvalues a = o"(e^) while the related eigenfunction has 
the form $(y) = f{my), i.e. n must be equated to in ([31), ([ll) and As an example we examine 
the profile W{s) = sins + sin 3s + cos 2s. We calculate the concrete asymptotic expansion of the 
unstable eigenvalues and compare it with the numeric results. The comparison exhibits very good 
coincidence. 
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It should be noted that function 7(o"o) has been introduced originally in [3] in order to formulate 
the necessary condition for the long-wave instability of the channel flows. However, the result of 
[3] exploits the monotonicity of the shear flow profile that makes it inapplicable to the oscillating 
flows. At the same time, the asymptotic we found can be considered as long-wave one since the 
problem ([5]) in fact depends on cP'e^ only when n = 0. 

Dispersion equation for an abstract spectral problem. 

Let K = -ff (cr, e) be linear operator defined for every a in some neighbourhood of do and for every 
e in some neighbourhood of zero. Assume that domains both of H and of H* do not depend on 
a and e, the resolvent of K is compact for every a and e and depends continuosly on a and e in 
uniform topology. Consider the spectral problem 

//(c7,e)/ = 0, (11) 

where a is considered as spectral parameter. Assume that //(ctq, 0)*== -ffoo has non-trivial kernel: 
keri^oo 7^ {0}, i- e. do is an eigenvalue for e = 0. Let us reduce the spectral problem ([TT]) to a 
scalar equation for parameters a and e. 

Operator //qo has compact resolvent, therefore its kernel is finite-dimensional and dimker //qo = 
dimker f^QQ. Let A'^ = dimker Hqq^ Pi, P2- projectors (may be non self-adjoint) onto A^-dimensional 
subspaces ker ifgo and ker Hqq. Let Qk = I — Pk, k = 1,2. Set 

f = Pif + Qif = C + 0. (12) 

Apply P2 and Q2 to equation (fTTI) . Write 

P2HPi^ + P2HQie = 0, (13) 

Q2HPi^ + Q2HQie = 0. (14) 

Operator Q2H00Q1 bijectively maps ImQi n V{H) on Im(52- Therefore the inverse operator 
{Q2H{a,e)Qi)^^ : Im(52 ^ ImQi is bounded and continuous in a small neighborhood of (cJoiO). 
We resolve equation (I14p in 6 and then substitute the solution into the (I13p . This gives an equation 
in a finite-dimensional subspace ker Hqq 

P2Hi - P2HQi{Q2HQi)-^Q2Hi = 0. 

We equate the determinant of this equation to zero and arrive to a scalar equation. This is the 
dispersion equation for a and e. The structure of the asymptotic of (7(e) can be deduced from 
Newton diagram (see j^). 

The limit operator kernel 

The operator Hqq = H{ao,0) of the problem ([H) is defined by 

^00/ = {(70 + iaW)f" - iaW'f = 

ds \ ds a + laW J 
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For all (To it's kernel is non-trivial. It contains functions 

f{s) = Ciao + iaW) (15) 

with arbitrary constant C. If equation 



doesn't have solutions different from p5]) then one can apply the results jl]. We focus ourselves 
on the case of 2-dimensional kernel of Hqq i.e. we assume that equation (I16p has two linearly 
independent solutions. 

Integrating equation ([TB]) . we obtain 

d f Ci 



ds a + iaW (a + iaWy ' 
The solvability condition for this equation is 

Ci < (a + iaW)-'^ >= 0, 

where < .. > denotes average: < / >= ^ f^^ /(s) ds. Thus, the kernel of Hqq is 2-dimensional if 
and only if 

7(fTo) =< {a + iaW)-^ >= 0. (17) 
Then the general solution of (jl6p has the form 



(fiis) = (70 + iaW{s), <p2is) = (fi 

dy 



(fjo + iaW{s)) 



1 

° (18) 



{ao + iaW{y)f 



We note that 7 (do) is analytical function provided RecJo 7^ 0. Its zeros are the multiple eigenvalues 
we are looking for. 

The kernel of the adjoint operator 

Let ctq satisfy p!7|) . Consider the adjoint operator 

1 <i ,_.,n2* 



Let us integrate equation 

l(w-,Mr)^|).„, 



and then divide the result by (ctq — iaW)^, and then integrate once more. As a result we have 

dy 



s 

/ = Ci + C2 




{a*Q-iaW{y)) 



2 ■ 
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According to (jl7p . / is periodic for all Ci and C2. Therefore, functions 

s 

dy 



Ms) 



^^^^^ =/k 



iaW{y)f 



(19) 



span the kernel of HZ 



00- 



Dispersion equation for the problem 

Let us construct dispersion equation for the case of 2-dimensional kernel of limit operator. Let Pi 
and P2 be the orthogonal projectors onto subspaces keri^oo ^-nd keri^QQ . Functions V'l and ^2 are 
orthogonal so that 

P2j = -ri — TT^i + Ti — rT^2 =</> +7- — 

Although and (p2 are not orthogonal, it is clear that projection Pi/ = if and only if (/, (pi) = 
and (/, Lp2) = 0. We define projectors (5i,2 setting (5i,2 = I — Pi,2- 
Let = Qif. Then into sum 

d = piO^^^ + (320^^\ (20) 
where functions ^'^-'^'^^ satisfy an equation 



Q2HQie^^ + Q2H^i = 0, i = l,2. 



(21) 



(this representation can be seen from (jl3p ). Then the equation (jl4p in 2-dimensional subspace 
ker Hqq is equivalent to 

/?l(/7((^i + e«),V'l) +/?2(^(9^2 + e^'^),V'l) =0, 
Pl{H{ipi + e'^^^),^2)+P2{H{ip2 + 9^^^),^2) =0. 

The equating of the determinant of this system to zero gives us the dispersion equation 



A = 



(p'(vPi + eW),Vi) {H{ip2 + e(^^),i>i) 

(P'((^i+e«),^/>2) {H{^2 + 0^^^),iJ2) 



0. 



(22) 



The construction of the asymptotic 



Let's plot Newton diagram to determine the structure of the eigenvalue asymptotic (see [3]). We 
expand function 9^^''^\ operator H and determinant A in the powers of e: 



0» 



H2f 



H{a,e) = Ho{a) + e^Hi{a)+e'^H2{a) 

2fr ' ' 



'f- 



; Hif = 2in(Pf'; 



(23) 
(24) 

(25) 



-(a^ + n^)^/; ^{s) = a + iaWis). 
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A = Ao + e^Ai +e2A2 + ... 
According to expansion (I24|) equations ()2ip take the form 

= -Q2Hoipi - e^Q2HiLpi - e^Q2H2^Pi. 

The substitution of ([23]) gives series of equations for coefficients 9^l\ i = 1,2 

Q2Ho9q'' = —Q2Hoipi, 
Q2Ho6f + Q2Hief = -Q2Hiipi, 

Q2Hq92^ + Q2Hi9i^ + Q2H29q^ = —Q2H2^i, 

Q2HQef + Q2Hiefl^ + Q2H2ef_2 = o. 



(26) 



k = l 
k = 2 
k>3 



(27) 



Theorem 1. a) Coefficients Ao(o") and Ai{a) identically equal to zero. 

b) If n ^ then A2((To) 7^ 0, and such (Tq produce no analytical branches of eigenvalues cr(e). 

c) lfn = then A2(fTo) = and A2fc+i(<T) = 0, /c G N. 

d) Letn = and J^^ ds / 0. // (To is real then assume additionally that ^ < >. 
Then ao is a limit point for some branch of simple eigenvalues cr^e"^), and 



where A4(fTo) and ^^(o"o) are defined by equalities 



0-2 



A4(ao) 

dA2 I „ \ ' 



(28) 



dA2 



da 



[H2(pi,1pl} {H2ip2,'4'l) 
[H2<^l,1p2) (H2ip2,1p2) 



2o? <ip\><ip{^ > . 



To prove the theorem we have to expand the scalar products (yHnpj, tpk) and {^HiOm ,'^k) in the 
powers of e. 

Lemma 1. For all a: Reo" 7^ following statements hold: 

a) {Hofj,iJi) = {Ho9^ri\4'i)=0; 

b) {Hnpi,ipi) = 0; c) {Hiip2,'4>i) = in; 
d) {Hiipi,tp2) = -in; e) (i/o¥'i,'02) =0 

fj If a ^ ao then any solution o/(52-f^o^ = can be represented in the form = Ci4) + C2S,, where 
^ - particular solution of HqS, = V'2) function (f) defined in (|25p . 

Proof of the lemma. Equahties a,b,c,d and e are straightforward. To get statement f one 
can use the equivalence of equations Q2H0C = and Hq^ = aiip2 + a.2'4'2- Then ai = (according 
to the solvability condition). Then the general solution of homogeneous equation Hq^^ = is Ci(/) 
while Q2C is particular solution of the equation Hq£^ = 02^2- This completes the proof. 

Proof of the theorem 1. 
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a) According to statement a of lemma 1 
Ao 



+ 



+ 



+ 



{H,{ipi + 9l^^),iPi) (ifi(^2 + ef ),^i) 
{Hoi^i + ^ ^2) {Ho{v2 + ), ^2) 

(iJo(<^l+e?^^2) (iJo(<^2 + ^S'^),V'2) 

(i7o(v^l + ^S'^ V'l) (i^0('^2 + ), V'l) 

{Hii^i + ei^^),i;2) {Hi{y^2 + ), V'2) 



+ 



+ 



+ 



(2)n 



{Hoivi + e'i^),^2) {Ho{^2 + ), V'2) 
According to statement f of the lemma + ^J-) = C}"*)0 + Cf)e Then 

{H,{^m + V'2) =ct^ {Hoi i^2) =ct^ (V2, V2), 

(the statement b of lemma 1 was used for the last equality). Thus, we get a determinant with 
proportional rows 

Ci'\i^2,^2) (V'2,V'2) 

b) Writing coefficient A2 similarly to Ai we obtain the sum of eight determinants. Half of them 
are equal to zero according to the statement a of lemma 1. As a result we arrive at equality 



Ai = 



= 0. 



A, 



+ 



+ 



+ 



(i/2(¥'l + tf^V'l) {H2{<P2+0^o^),4^l) 

(/f0(¥'l + ^S'^^2) (/^o(¥'2+ef ),^2) 

{H,{<fi + ^S'^ V'l) {Hi{<P2 + ), ^1) 

{H,{<fi + ^S'^^2) {Hl{<f2 + ), ^^2) 

{Ho9'^\^P2) {Hoef\^2) 



+ 



+ 



(29) 



+ 
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Let a = (To. Then the first, second and forth determinants are equal to zero as Imi^oo -L V'l 2- since 



-f^oo</'i,2 = equations (f27l) yield ' 
lemma 1 we get 

A2(ao) = 

4k 



-n 



4 
3 

2# 

1 




a) 



provided a = a^. Using statements b,c and d of the 

in 
-in {Hnp2,tp2) 

b) 



m 




m 



Fig. 1. Newton diagram 

If n 7^ then A2(cto) ^ 0. The related Newton diagram is presented on fig. la. Line k = 
corresponds to the coefficient Aq, line k = 1 — to the coefficient Ai and so on. Since A2(o"o) 7^ 0, 
Newton diagram includes point (0,2), while there are no points below the line k = 2 (due to 
statement a of theorem 1). Therefore, equation (j22p doesn't have solutions nearby do when e — > 0. 

c) If n = then Hi = and operator H in fact depends on e^. Consequently, coefficients A2a,.+i 
are equal to zero. 



(1 2) 

Let's calculate A/i{(Tq). Taking into account that 6q ' (uq) = and Imffoo -L ^"1,25 we get 

[H2ipi,tpl) {H2(p2,i'l) 
{H2ipi,1p2) (-f^2952, V'2) 



A4(ao) 



Let's calculate ^^(cjo). Taking into account that Hi = and (|29|) we get 



dA2 

da 



{H2ipi,tpl) {H2(p2,1pl) 



'dHo 



Then (^^</'i,'02) = by statement e of the lemma. Since = f", we have 



-^i^o) = (^2'^l,V'l)('/2'2,V'2). 



(30) 



Let's determine scalar products (ff2V'iiV'i) aiid (932, ■02). 



{H2ipi,i^i) = -— J ^p\ds = -Q^[Re((T^ 



re ^im I 



If (Tq.^ 7^ 0, then imaginary part of the scalar product is nonzero otherwise we have to require that 
(Tq 7^ < W"^ > in order ensure that ^^(cq) is nonzero. Let's transform the second multiplier in 
formula (1301) 



¥'2>'02)=^ \^2[^] ds = - lipi( /--2 



IT 



1 



ds- 



ds 



iff 



ds. 
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Assuming that < ip^'^ we arrive at the Newton diagram presented on fig. lb which, in turn, 
imphes (I28p . The proof is completed. 



Computing experiment 

Consider flow profile W{s) = sins + sin 3s + cos 2s. If n = 0, then problem ([4]) depends on 
only one parameter: after substitution a = iaa, e = e/a only one parameter e and unknown 
eigenvalue a left. Therefore, we can set a = 1 with no losses in generality. Function 7 defined in 
([T7D vanishes in the point ctq = 0.3543 — 0.6366i. Coefficient (T2 of asymptotic expansion ([25]) is 
(72 = —0.2568 + 0.2393i. The asymptotic aap = cfq + a2S^ + ■ ■ ■ demonstrates very good agreement 
with the Numerical results (see fig. 2 and 3). 
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Fig. 2. Plots of the numerical (solid line) and asymptotic (dashed line) solution (t{£). Abscissa 
is e, ordinate axises are Reo" and Imo" respectively. 



0.15 



0.1 



0.05 



\cr — (j 



ap I 



0.2 



0.4 



0.6 



0.8 



8 



Fig. 3. The relative difference between the numerical solution and asymptotic. Abscissa is e, 

ordinate axis is 
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